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ABSTRACT 
We explore a technique for identifying the highest redshift {z > 4) sources in Herschel/ SPIKE 
and BLAST submihimeter surveys by locahzing the position of the far-infrared dust peak. Just as 
Spitzer/IRAC was used to identify stehar 'bump' sources, the far-IR peak is also a redshift indicator; 
although, the latter also depends on the average dust temperature. We demonstrate the wide range 
of allowable redshifts for a reasonable range of dust temperatures and show that it is impossible to 
constraint the redshift of individual objects using solely the position of the far-IR peak. By fitting 
spectral energy distribution models to simulated 7/ersc/ie^/SPIRE photometry we show the utility 
of radio and/or far-infrared data in breaking this degeneracy. With prior knowledge of the dust 
temperature distribution it is possible to obtain statistical samples of high redshift submillimeter 
galaxy candidates. We apply this technique to the BLAST survey of ECDFS to constrain the number 
of dusty galaxies at z > 4. We find 8 ± 2 galaxies with fiux density ratios of 5500 > S350; this sets an 
upper limit of 17 ± 4 deg~^ if we assume all are at z > 4. This is < 35 % of all 500 //m-selected galaxies 
down to S500 > 45mJy (Ljr > 2 x IO^^Lq for z > 4). Modeling with conventional temperature and 
redshift distributions estimates the percentage of these 500 /xm peak galaxies at z > 4 to be between 
10-85%. Our results are consistent with other estimates of the number density of very high redshift 
submillimeter galaxies and follows the decline in the star formation rate density at z > 4. 

Subject headings: galaxies: evolution — galaxies: starburst — galaxies: high-redshift — infrared: 
galaxies — submillimeter — techniques: photometric 



1. INTRODUCTION 

The prevalence of dusty, infrared luminous galaxies at 
early times tests our theories about early galaxy forma- 
tion and the production of metals. Semi-analytic mod- 
els predict very small numbers of submillimeter galax- 
ies (SMGs) at extreme redshifts (z > 4, e.g. Baugh et 
al. 2005), and the best observational constraints on the 
redshift distribution of SMGs are consistent with this 
(Chapman et al. 2005; Swinbank et al. 2009). Much 
attention has been given to recent studies which have 
spectroscopically confirmed a small number of SMGs at 
z > 4 (Capak et al. 2008; Coppin et al. 2009; Daddi et 
al. 2009ab; Knudsen et al. 2010); however, these sources 
are still only a small fraction of the submillimeter pop- 
ulation. Small number statistics and cosmic variance in 
current submillimeter surveys prohibit an accurate mea- 
sure of the space density of the most distant SMGs. 

The question of the prominence of the high redshift tail 
of dusty SMGs can be addressed with imminent submil- 
limeter wide area surveys such as those with the Herschel 
Space Observatory (Pilbratt 2001), the SCUBA-2 cam- 
era on the JCMT (Holland et al. 2006), and the AzTEC 
camera on the Large Millimeter Telescope (LMT, Wil- 
son et al. 2008). The Herschel Space Observatory is 
currently collecting data from 70-500 /im over wide sur- 
vey fields. With an aperture of 3.5 m, Herschel is only 
sensitive to ultra-luminous infrared galaxies (ULIRGs, 
L/fl > 10'^^Lq) above z > 2 at 250-500 //m before it 
reaches the confusion limit. One advantage of SPIRE 
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surveys is the simultaneous observations at three wave- 
lengths; 250, 350 and 500 /im, which samples the peak of 
the far-IR spectral energy distribution (SED) due to dust 
emission in high redshift galaxies. With prior knowledge 
of the dust properties which describe the shape of the 
SED, the position of the far-IR peak is sensitive to the 
redshift of the galaxy. 

In this paper, we discuss the power and limitations of 
using the far-IR peak as a redshift indicator in submil- 
limeter surveys. We apply this technique to the current 
BLAST survey of ECDFS and place an upper limit on 
the number density of 500 /xm-selected SMGs at z > 4. 

Throughout this paper we assume a standard cos- 



mology with Hq 
Qa = 0.73. 
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2. FAR-IR 'BUMP' TECHNIQUE 



The 1.6 /im bump in the near-infrared SED of galaxies 
arises due to a minimum in the H~ opacity in the spectra 
of cool stars. This bump has been used in extragalac- 
tic Spitzer/IRAC surveys to identify sources in a specific 
redshift range (e.g. Wright et al. 1994, Sawicki 2002, Far- 
rah et al. 2008). A similar, but factor of '^ 2 — 3 broader, 
peak exists in the SED of galaxies at far-IR wavelengths 
which is caused by the integrated thermal emission from 
dust of different temperatures. This SED is often pa- 
rameterized as a blackbody distribution of a particular 
far- infrared color temperature (T^ust) with an additional 
term to account for the dust emissivitjQ (e.g. Blain et 
al. 2002). 

The challenge with using the far-IR bump as a redshift 
indicator is that the average dust temperature and the 

* In this paper, we assume a dust emissivity, /3 = 1.5. The choice 
of emissivity affects the derived color temperature. 
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redshift are degenerate (e.g. Blain 1999; Blain, Barnard 
& Chapman 2003). Wien's displacement law tells us that 
the wavelength of the peak of a blackbody scales with 
the dust temperature. Coupling this with the displace- 
ment of the peak due to redshift we find that the ob- 
served wavelength of the far-IR peak depends linearly on 
both the redshift and the inverse of the dust temperature: 
AX/ oc (1 + z)/Tdust. 

The most luminous galaxies in the local Universe 
show far-infrared color temperatures that peak at '^40 K 
although there is a significant scatter in luminosity- 
temperature parameter space (e.g. Dunne et al. 2000). 
For a reasonable range of dust temperatures for ULIRGs, 
sources that peak at 500 /xm can be found anywhere be- 
tween z '^ 3-6. This is a very large range in redshift and 
without prior information about the dust temperature 
it is difficult to further constrain the redshift with the 
SPIRE data alone. 

3. HIGH REDSHIFT CANDIDATES IN THE ECDFS BLAST 

SURVEY 

In order to test this technique on observations we use 
the data from the BLAST survey of ECDFS (Devlin et 
al. 2009). BLAST has the same three detectors as Her- 
sc/ie//SPIRE but with beam sizes that are twice as big; 
36, 42, 60 arcsecs FWHM at 250, 350 and 500 /xm respec- 
tively (Marsden et al. 2009). We use the publicly released 
BLAST mapf0 and the matched filter catalogs (Chapin 
et al. 2010). These matched filter catalogs do a much 
better job at de-blending adjacent sources and result in 
higher signal-to-noise ratios than previous BLAST cata- 
logs (see Chapin et al. 2010 for further details). We re- 
strict our analysis to the central 0.47 deg^ of the BLAST 
image where crinst < 10 mJy at 500 /im. We have con- 
servatively added the confusion noise in quadrature with 
the instrument noise (cTconf + ^inst = '''tot) foi" all anal- 
ysis in this paper. We only consider detections if they 
are > 5crinst (> 3crtot)- In this central region there are 23 
galaxies robustly detected at 500 [iva. with 6*500 > 45 mJy. 
We match the 500 /im-selected galaxies to the matched 
filter catalogs at 250 /im and 350 /im; we consider any 
> 5(Tinst detection within a radius of 60 arcsecs to be the 
counterpart to the 500 /im emission. For sources which 
are undetected at these other wavelengths we set upper 
limits on the 350 and 250 /im fluxes from the 90% com- 
pleteness limits of the survey (Chapin et al. 2010) which 
is equivalent to 3crtot- 

Of these 23 500 //m-selected galaxies, 8 have 5*500 > 
5*350 indicating that the dust SED peaks at or near 
500 /im - we refer to these galaxies as 500 /im 'peak- 
ers' (although their peak may be at even longer wave- 
lengths). Interestingly, none of these 500 /im peakers (de- 
tected at > 5(Tinst) are detected above 5crinst at 250 and 
350 /im. We inspected all of these candidates by hand 
in the BLAST images to make sure they are not obvi- 
ously blended. We have not attempted to correct for 
flux boosting (e.g. Coppin et al. 2006) in this analysis 
but, since we are not interested in the absolute fluxes 
and only the relative colors, this is less of a concern. In 
addition by considering only higher signal-to-noise ratio 
sources we ensure that this effect is minimal. In order to 
assess how photometric uncertainties (from both confu- 

^ http://blastexperiment.info/rcsults.php 



sion and instrument noise) affect the number of 500 /im 
peakers identified, we ran a Monte Carlo simulation to 
randomly sample the fluxes of each of the 23 sources 
from a Gaussian distribution S ± Ctot ■ We find that pho- 
tometric uncertainties can vary the number of sources 
with ^500 > ^350 by 2. Thus our final number of 500 /im 
peakers is 8 ± 2. 

These 8 high redshift candidates are listed in Table 1 
along with their fiuxes and Fig. [1] shows an example SED 
for one source. Five of these 500 /im peakers are within 
the ECDFS field and one is within the smaller GOODS- 
S region (see Table 1). The source within GOODS-S is 
associated (6 arcsecs away) from an AzTEC 1.1mm de- 
tected source (GSll, Scott et al. 2010); however, fiux 
boosting in the BLAST bands prohibit a detailed com- 
parison. 

Knowing that a source peaks at (or near) 500 /im does 
not provide a definitive redshift due to the degeneracy 
between redshift and temperature (e.g. Fig. [5]). If we 
conservatively assume that all 500 /im peakers are at z > 
4 we can obtain an upper limit on the number density 
of z > 4 galaxies predicted in Herschel SPIRE surveys. 
Based on our analysis of the BLAST data, the number 
density of galaxies (with S500 > 45 mJy) which peak at 
500 /im is 17 ± 4deg~^. We expect that some of these 
500 /im peakers will be at z < 4 since Fig. [2] shows that 
sources cooler than 40 K can have 6*500 > 5*350 but be at 
z < 4. 

To get an estimate of the contamination of low red- 
shift, cooler galaxies in the 500 /im peaker sample we 
must assume a redshift and dust temperature distribu- 
tion. Assuming a flat redshift distributioip and a dust 
color temperature distribution of 35 ± 7 Klj, we estimate 
that 15 % of 500 /im peakers will be z < 4 (left panel of 
Fig. [2). In addition, we must also account for the con- 
tribution from warmer sources at z > 4 but which peak 
at lower wavelengths (6500 < S350). Assuming a sam- 
ple of sources at z > 4 again following a dust tempera- 
ture distribution of 35 ± 7K, we find that only 10 % have 
5500 < 5350. These two factors roughly cancel each other 
and we are left with a constraint on the number density 
of z > 4 sources (with S500 > 45mJy) of < 17deg~^. 

The above limit on the number density of z > 4 SMGs 
ambitiously assumes a flat redshift distribution. If in- 
stead we assume the observed redshift distribution of 
SMGs (2.2 ±0.8, Chapman et al. 2005; Pope et al. 2006), 
we find only 10% of the 500 /im peakers at z > 4 (right 
panel of Fig. [2]) which brings the number density esti- 
mate down to 2deg~^. Based on our analysis of the 
BLAST data, we conclude that the number density of 
S500 > 45mJy sources at z > 4 is < 17dcg^^ and could 
be as low as 2 deg^^. 

With some assumptions about the dust temperature 
distribution, we have shown that the submillimeter color 
selection, 6500 > 6350 is a plausible way to identify sam- 
ples containing the highest redshift SMGs selected at 
these wavelengths (albeit these samples contain contam- 

® A flat redshift distribution has an equal number of objects per 
redshift interval (constant dn/dz). This is roughly equivalent to a 
non-evolving luminosity function at 2 > 2. 

^ This dust temperature distribution is consistent with local 



ULIRGs (Dunne et al. 2000) as well as z 
et al. 2005; Pope et al. 2006). 



2 SMGs (Chapman 
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ination from lower redshift objects). Without further 
multi-wavelength data, it is impossible to tell individu- 
ally which of the 8 candidates are at z > 4, and which 
are just peaking at 500 /im because they have cooler dust 
temperatures. One way to test the reality of these high 
redshift candidates is to further constrain their SEDs 
with data at longer submm wavelengths such as those 
from the 870 ^m LABOCA survey of ECDFS (Weiss et 
al. 2009). In order to do a robust comparison between 
the LABOCA and BLAST fluxes, a correction for flux 
boosting needs to be applied to both surveys using the 
same method (e.g. Coppin et al. 2006). This is beyond 
the scope of this paper but should be possible with full 
access to the signal and noise maps from both surveys. 
Another potential way to break the degeneracy between 
redshift and temperature for these candidates is with the 
radio and/or far-IR (e.g. 100 /xm) flux; we discuss this 
further in the next section. 

4. USING RADIO AND FAR-IR DATA TO BREAK THE Z-T 
DEGENERACY 

Photometry from 250-500 /im alone can provide an ac- 
curate measure of (l-|-z)/Tdust- Fig.[3]shows a simulation 
of Herschel SPIRE photometry for a source that peaks at 
500 /zm. Fitting modified blackbody models to this data 
we find the pairs of T and z that fit these data (panel b 
of Fig. [3]) and the resulting LpiR (panel c of Fig. [3|). 

The radio emission from galaxies is known to corre- 
late well with the infrared emission in the local Universe 
(Condon 1992). Because of this radio-infrared correla- 
tion, radio data is often used to help identify counter- 
parts to SMGs (e.g. Ivison et al. 2002), and to constrain 
the redshift (e.g. Carilli & Yun 1999; Hughes et al. 2002; 
Aretxaga et al. 2007). We demonstrate the latter point 
in panel d of Fig. [3l Here we assume the average value 
and scatter in the q parameter from fitting local galaxies 
(Yun, Reddy & Condon 2001). The radio flux is not de- 
generate with the redshift; however, unfortunately, the 
large scatter in the local radio-IR correlation provides 
fairly loose constraints on the redshift of individual ob- 
jects from the radio flux. Nevertheless, a large fraction of 
the bright SMGs detected with SPIRE will be above the 
flux limits of deep (cTradio < 5 /J.Jy) 1.4 GHz radio surveys. 
With large statistical samples, the redshift distribution 
inferred from including radio observations should be ro- 
bust, assuming the radio-IR correlation holds out to high 
redshift. 

Below ^ 50 fiui rest-frame, the SED of galaxies is no 
longer dominated by the single temperature modified 
blackbody but instead by warmer dust components. This 
can be seen in panel a of Fig.[3]where the dash-dot curve 
is a Chary & Elbaz (CEOl) template representative of 
local star forming galaxies. Fitting the CEOl templates 
(allowing the luminosity and temperature to vary) to the 
simulated SPIRE photometry we estimate the 100 /im 
flux as a function of redshift (panel e of Fig. [3|). While 
the curve is not smooth (due to discrete templates within 
the CEOl library), overall the far-IR flux can constrain 
the redshift for this simulated galaxy with no prior as- 
sumption on its dust temperature. The dotted line in 
panel e of Fig.[3]shows the sensitivity of the deepest Her- 
schel PACS extragalactic survey (GOODS Herschel, PI: 
D. Elbaz); this deep survey should detect bright SMGs 
beyond z ~ 3 by measuring the redshifted emission from 



warm dust that is heated in O&B star-forming regions 
and appears to correlate well with ongoing star-formation 
(Calzetti et al. 2005). 

5. DISCUSSION 

From the BLAST ECDFS data we estimate a space 
density of < 17deg^^ for sources with S500 > 45mJy 
at 2 > 4; this is < 35 % (8/23) of all sources selected 
down to these depths in the BLAST data. These bright 
500 /xm-selected SMGs at z > 4 contribute < 1% of the 
total background emission at 500 ^m (Fixsen et al. 1998). 

Constrained by a recent analysis of the infrared ex- 
tragalactic background light (EBL), galaxies evolution 
modelspredict the source density of galaxies with S500 > 
27mJj|3at z > 4 to be 1-5 deg'^ (Chary & Pope 2010). 
The range of predicted number density comes from con- 
sidering two evolutions of the models both of which fit 
the observed EBL. We find that our current observational 
constraint is consistent with both evolutionary scenarios. 
Future submm surveys with much larger telescopes such 
as JCMT, LMT and CCAT wiU reach much deeper flux 
limits before hitting the confusion limit. With these sur- 
veys, we can expect to detect many more dusty galaxies 
at z > 4; the models in Chary & Pope (2010) predict 
a number density of 40-200 deg~^ at z > 4 for galaxies 
down to 5 mjy at 500 /im. 

In the past two years, the tally of z > 4 spectroscopi- 
cally conflrmed SMGs has grown from zero to flve0 (Ca- 
pak et al. 2008; Coppin et al. 2009; Daddi et al. 2009a,b). 
Coppin et al. (2009) estimate the number density of 
SMGs at z > 4 by combining these 5 sources from 3 inde- 
pendent submillimeter surveys and find a lower limit of 
> 7deg~^. While this is consistent with our predictions 
for 500 /^m-selected galaxies we note that the sensitivi- 
ties of the BLAST and submillimeter surveys are not the 
same at z; = 4; only 1/5 of the z > 4 spectroscopically 
confirmed SMGs might be detected down to the BLAST 
ECDFS survey depths (assuming a typical SMG dust 
temperature). Folding this in, the Coppin et al. (2009) 
estimate down to Ssso > 20mJy scales to > 1.4deg~^ at 
z > 4. 

There is some evidence that large numbers of very mas- 
sive galaxies at z > 4 pose a substantial challenge to stan- 
dard galaxy formation theory (e.g. Baugh et al. 2005). 
Fig. [4] summarizes the current observational constraints 
on the number density of bright SMGs as a function of 
redshift. Blue and green bars are for 850 /xm-selected 
SMGs using values from the literature (Chapman et 
al. 2005; Coppin et al. 2009; Pope et al. 2006) and the 
red bar is our new upper limit from the analysis in this 
paper of the BLAST 500 /Ltm-selected sources. Overall, 
there does not appear to be room for a large fraction of 
the SMG population to be above z = 4 down to these 
flux limits. 

6. SUMMARY 



* In order to compare to the models we need to assume a rough 
correction for flux boosting of the BLAST measurements. Based 
on Eales et al. (2009), we adopt an average correction of 0.6 to 
the BLAST fluxes which effectively brings the flux limit down to 
27mJy. 

^ A sixth 2 > 4 SMG recently appeared in the literature (Knud- 
sen et al. 2010). 
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We have shown that the degeneracy between dust tem- 
perature and redshift severely hmits the use of the far- 
IR dust peak alone as a redshift indicator. The color 
selection of 6*500 > 'S'350 (and 5*500 > '5'25o) can be used 
to identify samples of candidate high redshift galaxies. 
In order to further constrain the redshift of individual 
sources requires prior knowledge of the dust tempera- 
ture or additional multi-wavelength data, particularly in 
the radio or at 100 ^m. 

We present 8 candidate 500 /im 'peakers' from the 
BLAST ECDFS survey with S500 > 45 mJy; a fraction 
of which we expect to be at z > 4 depending on the dis- 
tribution of dust temperatures. The number density of 
these high redshift candidates is < 17deg~^ and is con- 
sistent with the number densities of the brightest 850 /im- 
selected galaxies. The corresponding fraction of 500 /ini- 
selected galaxies in the BLAST survey which could be at 



2 > 4 is < 35%. 

In order to further constrain the space density of the 
most distant SMGs will require deep, wide area surveys 
with Herschel/SPIRE in addition to sufficiently deep 
multi- wavelength data to weed out the genuine high red- 
shift candidates from 500 /im 'peakers' that are lower 
redshift galaxies with cooler than average dust tempera- 
tures. Regardless, our upper limits on the number den- 
sity of z > 4 luminous dusty galaxies suggests a strong 
decline in their number density from z ~ 2. 

We thank the referee for constructive comments which 
improved the presentation of this work. We thank the 
BLAST team for making their data public. AP acknowl- 
edges support provided by NASA through the Spitzer 
Space Telescope Fellowship Program, through a contract 
issued by the Jet Propulsion Laboratory, California In- 
stitute of Technology under a contract with NASA. 
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Fig. 1. — Example far-IR SED of one of our high redshift candidates from the BLAST ECDFS survey. The current data and limits are 
consistent with very cold 2 ~ 2 sources as well as 2 ~ 4 sources with average dust temperatures. Stronger limits at 250 and 350 /.im will be 
required to alleviate this degeneracy. 
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Fig. 2. — Criteria for selecting 500 /jm peakers and the degeneracy between dust temperature and redshift. The curves show the S500/S350 
flux density ratio of an SED (single temperature modified blackbody with /3 = 1.5) for a given T^ugt (see color legend) as a function of 
redshift. The horizontal dotted line shows the criteria used to define a source that peaks at 500 ^m; there is a range of redshift and dust 
temperature combinations that produces this ratio. The data points show a Monte Carlo simulation of 2000 galaxies following a dust 
temperature distribution of 35 + / — 7K, and a uniform redshift distribution (left panel) or z = 2.2 + / — 0.8 (right panel). The uniform 
redshift distribution results in the majority of 500 /im peakers being at 2: > 4, whereas the Gaussian redshift distribution puts very few 
sources above z = 4 and the majority of 500 fim peakers at z < 4. 
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Fig. 3. — Example of degeneracy in single temperature modified blackbody SED fitting, a) Simulated SPIRE photometry of a 500 /im 
peaker (S500 = 40mjy). Flux uncertainty includes both confusion {a-^anf ~ 6mjy) and instrument (fTinst ~ 1 mJy) noise (the shaded 
regions in b), c), d), and e) refiect the errors due to this measurement uncertainty). The dashed curve is the best-fit modified blackbody 
where temperature and redshift combinations can be anything along the line shown in b). We have excluded temperatures below 15K since 
these are unphysical due to heating by the CMB. c) The resulting Lpijj of the best-fit modified blackbody model as a function of redshift. 

d) The predicted radio fiux of the best-fit modified blackbody model as a function of redshift (solid curve). The dotted curves illustrate the 
uncertainty resulting from the scatter in the local radio-IR correlation. The dash-dot curve in a) is a CEOl template which shows how the 
true SEDs of galaxies diverge from a single temperature modified blackbody at wavelengths shorter than around 50 fim in the rest frame. 

e) Estimated 100 fim fiux for best-fit CEOl SED as a function of redshift. The dot-dash curve shows the 5cr depth of the deepest PACS 
extragalactic survey, GOODS Herschel (PI D. Elbaz). Deep radio and PACS observations can help break the degeneracy between redshift 
and dust temperature for the most luminous dusty galaxies. 
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Fig. 4. — Histogram showing the number density of submiUimeter-selected galaxies as a function of redshift. Our new upper limit from 
the 500 ^m BLAST analysis is shown as the open red bar. The blue and green bars are estimates for SCUBA 850 /.im-selected galaxies from 
the literature where filled bars denote spectroscopic samples {z < 4, Chapman et al. 2005; 2 > 4, Coppin et al. 2009), and open bars denote 
photometric redshift samples (Pope et al. 2006). *Note that the BLAST sample has not been corrected for flux boosting, completeness, or 
spurious sources. The combination of these effects will act to decrease the true flux limit of this sample. 
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TABLE 1 
Catalog of all candidate 500 ^m 'peaker; sources. IDs and fluxes come from Chapin ft al. (2010); the uncertainties listed 

HERE include BOTH INSTRUMENT AND CONFUSION NOISE (crtot)- THE CONFUSION NOISE IS ASSUMED TO BE 12, 13 AND 15mJy AT 500, 350 

AND 250 /iM, RESPECTIVELY (ChAPIN FT AL. 2010). UPPER LIMITS CORRESPOND TO THE FLUX WHERE THE CATALOGS ARE 90% COMPLETE 

WHICH IS EQUIVALENT TO ROUGHLY 3(Ttot ■ THE FlELD COLUMN INDICATES IF THE SOURCE IS WITHIN THE ECDFS (E) AND/OR THE 

GOODS-S (G) AREAS. 



ID 5500 5350 5250 Field 

(mjy) (mjy) (mjy) 

BLASTJ033311-275611 
BLAST J033129-275548 
BLAST J033215-275030 
BLAST J033254-273303 
BLAST J033346-274217 
BLAST J033256-280105 
BLAST J033159-273515 
BLAST J033317-274926 



56 ±14 


<45 


< 60 


E 


53 ±14 


<45 


<60 


E 


50 ±14 


<45 


<60 


EG 


50 ±14 


<45 


<60 




50 ±14 


<45 


< 60 




48 ±14 


<45 


< 60 


E 


47 ±14 


<45 


< 60 




46 ±14 


<45 


<60 


E 



